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We study the conductance properties of a straight two-dimensional electron waveguide with an 
s-like scatterer modeled by a single 5-function potential with a finite number of modes. Even such a 
simple system exhibits interesting resonance phenomena. These resonances are explained in terms 
of quasi-bound states both by using a direct solution of the Schrodinger equation and by studying 
the Green's function of the system. Using the Green's function we calculate the survival probability 
as well as the power absorption and show the influence of the quasi-bound states on these two 
quantities. 
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I. INTRODUCTION 

Two-dimensional electron systems have been studied 
extensively over the past few years, both because it be- 
came feasible to construct such systems for example at a 
GaAs/GaAlAs interface at low temperatures and because 
the conductance was shown to be directly related to the 
transmission properties of th&. system. This relation is 
known as Landauer's formula,Eru 



(1) 



where F denotes the conductance and T the full trans- 
mission function (spin degrees of freedom are neglected). 
Eq. (|l|) results in a quantized conductance for straight 
channels. If one adds a single attractive i5-function scat- 
terer to such a straight waveguide, the combined effect 
of the scatterer itself and the backscattering off the walls 
leads to interesting phenomena. This model was first 
suggestedpl^ Datta et al.a and later discussed by Bag- 
well et al.Q'El In the following we will take a closer look at 
the resonance phenomena this system produces. 



II. DIRECT SOLUTION OF THE SCHRODINGER 
EQUATION 

Let us first obtain a solution of the Schrodinger equa- 
tion for an electron in a two-dimensional waveguide with 
a (5-function scatterer. The Hamiltonian is given by 



H 



P 

2m 



+ Vix,y) + Vciy). 



(2) 



Vc represents a confinement potential restricting the 
movement of the electron to the range < y < DU The 
attractive scattering potential is given by 



V{x,y) = -fS{x)S{y ~yo) (7 < 0). 



(3) 



We can now expand any stationary solution ipEix^y) of 
the Schrodinger equation, HipE{x,y) = EipE{x,y), in 



a Fourier series with x-dependent expansion coefficients 
using the complete set of transversal modes, 



V'B {x,y) = '^Cn{x)xn{y), 



(4) 



n=l 



where Xn{y) — \J % sin ^y. Inserting this series into the 
Schrodinger equation and employing orthogonality of the 
transversal modes we obtain a set of coupled equations. 



dx^ 



where Af„ 



(^) 



y^Mm„Cn(a:)(5(a;), 



(5) 
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denote the cou- 
the wavevec- 



pling constants and kn — w ^2 jj2 
tor (Imfcn > 0). Away from the scatterer, x 0, the 
wavefunction must have the free-electron form. 



A pik„x I f> -ik„ 



X < 
X > 



(6) 



As -0 must be continuous at x = and its derivative must 
have a finite jump there, the same conditions must hold 
for the expansion coefficients Cn{x). Thus using these two 
conditions on Eq. (||) with the ansatz in Eq. (|^) yields 

A. + C„ + (7) 

= Mnm {Am + B^) . (8) 

m 

If ip is an evanescent mode, we can set kn — inn and 
must require An = and Z3„ = to have a normalizable 
wavefunction. The transmission coefficient for propagat- 
ing modes is then defined as Tmn 
total transmission function as 



t^gp and the 



(9) 
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where the sum extends over all propagating modes. The 
conductance is finally calculated using Eq. (|l|) . The set of 
Eqs. and (^) can only be solved numerically for a finite 
number of modes. We use the parameters D = 300 A for 
the width of the channel and yo = foi' the transver- 
sal position of the scatterer, the mass m — 0.067TOe as 
the effective mass of an electron in GaAs and a scatterer 
strength of 7 = — 7feVcm^. With these parameters the 
first four energy subbands (transversal modes) open up 
at El = 6.24 meV, E2 = 24.94 meV, E3 = 56.12 meV 
and Ei = 99.78 meV. Numerical results of Eq. (|l|) for a 
total number of modes rij = 6 and nt = 100 are shown 
in Fig. 1^. It shows the expected step-like behavior at 
every energy subband edge. Furthermore, just before a 
new subband opens up, it shows an interesting dip. The 
waveguide blocks transmission in the lower mode just be- 
fore a new higher mode opens up. It becomes completely 
opaque just before the second mode opens up. The drop 
in the conductance to the n-th level just before the next 
higher mode n + I opens up is in fact a cumulative ef- 
fect of all transmitting modes 1 through n as-jcan be seen 
from the individual transmission coefficients.l3 These dips 
correspond to a resonance structure of the system. They 
can be attributed to quasi-bound states of the system 
with a fwite lifetime, represented by a wavefunction of 
the fornJa 



V'(x,t) = tp{x)e 



ERt/ih Eft/h 



(10) 



where Eji and Ej denote the real and imaginary part 
respectively. These states are characterized by having 
a scattered wave even without an incident wave, i.e. 
l/7r„„ = 0. We thus look for corresponding poles of 
the transmission function near each of these resonances. 
Figures || (a) and (b) show the poles corresponding to 
the first two dips of Fig. |l| confirming our assertion. The 
poles are always located in the analytically continued sec- 
ond sheet of the square root function (ImA:„ < 0). For 
the plots in Fig. ^the square root function was chosen to 
have its branch cut on the negative imaginary axis which 
is visible as a discontinuity in the plots (It is thus the 
square root function with values on the first, "physical" 
sheet in the first through third quadrant of the complex 
energy plane and the square root function with values on 
the second sheet in the fourth quadrant that is used in 
Fig. I). 



III. GREEN'S FUNCTION APPROACH 

The resonances observed in Section || can be dete 
mined from the Green's function of the system as well 
We thus try to find the exact Green's function of the 
Hamiltonian H. This on the one hand allows for a direct 
and exact calculation of the quasi-bound state energies 
as poles of the Green's function (on the second sheet). 
On the other hand, the Green's function calculation is 
less compute-intensive and thus allows the calculation of 



both the survival probability and the power absorption 
of the waveguide (see Sec. IV and 0). We first obtain 
the solution of the Green's function equation for the free 

waveguide with Hq = + Vc{y)^ 



(11) 



^ikn \x — x' 
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kn is given as in Section ||. The retarded solution (-1-) 
is denoted by G^'^ , the advanced solution (— ) by G°"^. 
For the special potential V{x, y) the integral equation for 
the full Green's function can be solved to yieldMa 

G^'^ (x, y, x', y', E) = G"'«/^ (x, y, x' , y\ E) (12) 
0, E) G^-^/^ (0, yo, x^ z/^ E) 
l/^_GO^flM(o,2;„,0,zjo,i5) 

The transmission- is solely determined by the retarded 
Green's function.!!! We can thus read off the condition for 
the poles to be 



1 4m 



1 
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2/0 I TTT^ = 0. 



(13) 



It can be explicitly shown from Eq. (g_3|) that the poles are 
always located on the second sheet and hence the branch 
cut has to be chosen as described in Section ||. Evaluat- 
ing Eq. ([1^) again numerically for nt — 6 and nt = 100 
modes gives the pole locations in Table |. Apart from 
the quasi-bound states the attractive (S-function scatterer 
also exhibits one single bound state. The pole locations 
agree exactly with the results from Section |^ thus con- 
firming the interpretation of the resonance phenomena 
to quasi-bound states. For the discussed (5-function im- 
purity, the poles depend sensitively on the number of 
modes included in the computation. The typical behav- 
ior of the pole location vs. number of modes is shown in 
Fig. H for the pole just below the second subband edge. 
The plot demonstrates that the location of the poles and 
thus the conductance properties of the discussed system 
are not converging with an increasing number of modes. 
It can also be seen that they are systematically shifted 
to greater imaginary values, displayed in the broadening 
of the resonance. 

The sum in Eq. ( p^ ) is in fact logarithmically diverging 
as each term in the sum is positive and approximately 
proportional to l/n for large n. The obtained results 
thus have to be always discussed for a specific number of 
modes. Nevertheless, the numerical results are in quali- 
tative agreement with the behavixix of an s-like scatterer 
as discussed by Kunze and LenkJlj They thus still serve 
their purpose as a useful model if the (5-function scatterer 
with a finite number of modes is interpreted as an s-like 
scatterer with finite width D /nt in j/-direction instead of 
a true (5-function scatterer. 

Before proceeding let us review the resonances in the 
context of the last two sections. They are formed by all 
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modes of the system collectively. Both propagating and 
evanescent modes (and therefore inter-channel coupling) 
are necessary. The evanescent modes are needed to build 
up a bound or quasi-bound state, whereas the propa- 
gating modes probe that state and therefore display the 
structure of the resonance. This can be seen clearly in 
Fig. ^ which shows the position of the pole causing the 
resonance just below the threshold of the second propa- 
gating mode. Keeping only two modes, n = 1 which is 
propagating and n = 2 which is extended and evanescent 
just below the threshold of the n — 2 propagating mode, 
is necessary to create the pole and hence the resonance. 
This can be seen clearly in the inset of Fig.^ where we 
show the conductance with only the n = 1 propagat- 
ing mode and the conductance when both the n = 1 
and n — 2 modes are present. There is no resonance 
when only the n ~ 1 mode is present. Thus, excluding 
the evanescent modes leads to unphysical kinks in the 
conductance near the threshold. The main contribution 
to any given resonance stems from the evanescent mode 
that is about to become propagating, because it has an 
infinite decay range at the band edge. The dip as well 
as the position of the pole are only slightly modified if 
more than this evanescent mode are included, but neither 
the dip nor the pole exist without evanescent modes. The 
remaining evanescent modes only alter the width and po- 
sition of the resonances and hence play a minor role. It is 
interesting to note that when more than one scatterer is 
present, evanescent modes play an even more important 
role because they can localize coherently over more than 
one scatterer. The case of many scatterers has been dis- 
cussed for waveguides in Ref. [l4|-p7| and for closed quan- 



tum systems in Ref. 18 



IV. SURVIVAL PROBABILITY 

Up to now the Green's function only allowed for a more 
compact formulation. We now demonstrate that it will 
also give new insight on other quantities like the survival 
probability. The survival probability is defined as 



(14) 



where \ipi) is the initial state at time t = and |i/'(0) is 
the propagated state, which can be computed from 



'(/'(x, t) — lim 



^ I I dz 



' — oo+ze 

e-'t'G^(x,x',z)^,(x')l , ^^ix,y), (15) 



however, the Green's function does not have any poles 
and only the integrations along the cuts, performed on 
different sheets, remain. We now demonstrate this case 
explicitly, for which we write the initial state in the fol- 
lowing form ipi{x.) — X{x)Y{y). Both X{x) and Y{y) are 
localized and do not contain an explicit z-dependence. 
For the integration along the n-th cut we define 



i^„"(x,x',t) 



(16) 



Although the integration runs on two different sheets, it 
is only relevant for the n-th square root, i.e. for fc„. Hence 
only one term contributes in the sum and we obtain 



i^r(x,x',t) 



(17) 



2raiT 
fit ' 



e" X*n{v)Xn{y')- 



Putting everything together we find for the survival prob- 
ability of the ballistic wire 



P{t) 



1 



2raiT 
U 



E 



e R 



dyY{y)xn{y) 



dx /da:'e ""'^" X*{x)X{x') 



(18) 



In Fig. H we plot the survival probability for an initial 
state which is Gaussian-like localized in the x-direction 
and has a mode expansion for the transverse part, i.e. it 
couples equally strong to all channels. We distinguish 
between short, medium and long time behavior. |_.For 
the short time behavior we observe the Zeno effect,E3 i.e. 
an initially non-decaying behavior. In the medium time 
regime we observe oscillations of all contributing chan- 
nels in the decaying probability. The long time behavior 
is dominated by an 1/t behavior. 

For the case of a J-function impurity in the wire we pick 
up poles which contain terms like e^^'^i^st/?* g^j^j therefore 
justify Eq. (|To|). Experimentally they can be detected by 
means of a Fourier transform and thus provide a tool to 
probe the complex spectrum. In addition we get contri- 
butions from the cuts, which do not cancel out, but are 
not dominating either. 



V. POWER ABSORPTION 



where the Green's function enters explicitly into the cal- 
culation. The integration runs in the complex plane as 
shown in Fig. ^ (solid line). We can deform the contour 
and perform the integration as indicated by the dotted 
lines. In this case we pick up poles as well as integrations 
along our choice of branch cuts. For the ballistic case 



Another quantity of interest is the power absorption 
(P) of a quantum wire, which can be calculated from mi- 
croscopic theory. It is related to the AC conductance 
G{oj) via (P) — P^j,jgL^G(cj) and hence allows us to 
make a straight connection between quasi-bound states 
and dips in the conductance. 
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An electric AC field of amplitude Eq is applied to a 
region of length L, wbich is symmetric around the scat- 
terer. It can be showriE3 that in linear response the power 
absorption of the ballistic wire behaves for small w's like 



{P) = 




(19) 



with k'j^ = \J kp — (ij)^- The quasi-bound states have 
a profound influence on the power spectrum. Start- 
ing from the microscmic expression P{t) = / dr'E{r,t) ■ 
(j(r,t)) one can showO that the quasi-bound states give 
Lorentzian-like contributions to the power absorption in 
these systems. However it is not clear whether these are 
positive or negative. We calculate the power absorption 
from 



{P) 



X {hi 

— TlLO 

X (Ei + huj J^(r') El 



dr Idv' E{x)E{x') / dEi 



Mr) 



El + huj 



(20) 



The microscopic current elements are evaluated in a scat- 
tering state basis. Our numerical results are shown in 
Fig. ^. The applied field is of length L = 100 nm, the 
scattering states are normalized on L' = 4000 nm. We 
clearly reproduce the analytical result for small to, which 
holds surprisingly well even for larger values. For the 
case of a ^-function impurity the spectrum is still domi- 
nated by the ballistic background, which is not surprising 
because the impurity is strongly localized. The influence 
of the quasi-bound states can also be seen, however it is 
quite Siftaller than for the system discussed by Na and 
Reichl.til In contrast to their system we see dips rather 
than peaks. This could be due to the fact that we excite 
into a quasi-bound state from a continuum state, which is 
why the dip shifts as a function of the Fermi energy. Na 
and Reichl however do excite from one quasi-bound state 
into another so that the signs cancel each other. In our 
system these resonances are too small to be observed, 
which can be attributed to the fact that the perturba- 
tion causing the quasi-bound state is much smaller (a 
(5-function with a flnitejjiumber of modes) than the large 
cavity in their system.Ej On the other hand we can now 
make a clear and unambiguous connection between the 
dips in the DC conductance and the quasi-bound states. 



VI. CONCLUSION 

We have considered a simple model for an electron 
waveguide. The conductance was calculated via Lan- 
dauer's formula and shown to exhibit resonance phenom- 
ena. These resonances were attributed to quasi-bound 



states of the system by both looking at poles of the trans- 
mission coefficients and the Green's function. We have 
shown that the delta-function potential together with a 
finite number of modes models an s-like scatterer. Fur- 
thermore, we have demonstrated how these quasi-bound 
states influence the survival probability and the power 
absorption of the system. 
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FIG. 1. Total conductance for a waveguide with a 
5- function scatterer of strength 7 = — 7feVcm'^. The con- 
ductance goes to zero just before the second subband opens. 





FIG. 2. Transmission poles for 7 = — 7feVcm'^. The poles 
just below the second (a) and third (b) subband edge are 
shown together with the corresponding zeros which appear as 
dips in these logarithmic plots. The zero in (b) is off the real 
axis as the transmission decreases to 1 instead of at this 
resonance (cf. Fig. |lj). 
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FIG. 3. Poles of G in the complex energy plane depend- 
ing on the number of total modes nt- The pole for rit = 100 
is located at i? = 24.888 — 0.615imeV and corresponds to 
the one shown in Fig. |^ (a). Inset: Conductance for nt = 2 
(solid) and nt — number of propagating modes (dashed). For 
E > E2 — 4Ei both are identical. 



FIG. 4. Integration path for the survival probability be- 
fore (solid) and after deforming the contour (dotted) in the 
complex a-plane. 
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FIG. 5. The survival probability P{t) as a function of time. 
The inset illustrates the Zeno effect for short times. 
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FIG. 6. Power absorption as a function of ui in units of 
2^jj 2 for the ballistic (dashed) and the impurity wire (solid). 

The Fermi energy is close to the quasi-bound state energy. 



TABLE I. Poles of G for a scatterer strength of 
7 = — 7feVcm and a total number of modes nt — Q and 
nt = 100 



Pole 


E [meV] (nt = 6) 


E [meV] {nt = 100) 





4.676 


-8.656 


1 


24.861 - O.OSli 


24.888 - 0.615i 


2 


55.804 - 0.144i 


55.117 - 1.96H 


3 


99.205 - 0.262i 


98.073 - 2.561i 
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